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BEAM  GENERATED  VORTICITY  AND  CONVECTIVE  CHANNEL  MIXING 


I.  INTRODUCTION 

Recent  theory  and  numerical  computation  at  NRL  have  resulted  in  a  theoretical  formalism  for 
estimating  persistent  flows  which  mix  cold  ambient  air  into  a  beam-heated  channel.  In  this  paper  we 
also  give  a  quantitative  calibration  of  that  formalism  using  detailed  two-dimensional  simulations.  The 
importance  of  determining  the  mechanisms  and  scaling  laws  governing  the  dynamics  of  a  gaseous  chan¬ 
nel  heated  by  a  sequence  of  pulses  is  clear.  One  important  application  is  the  fixation  of  atmospheric 
nitrogen  by  lightning. 

The  probable  existence  of  such  a  cooling  mechanism  is  well  documented  in  NRL  laser  and 
discharge  data  showing  heated  channels  breaking  up  in  a  number  of  circumstances. 1  Data  also  exists  on 
"turbulent*  cooling  of  multipulse  channels  from  Astron2  (Fig.  1).  The  "classical  theory"  curve  in  that 
figure  shows  that  the  thermal  conduction  cooling  rate  is  significantly  slower  than  observed  in 
experiments— pointing  to  a  convective  "turbulent"  mechanism.  Any  reasonable  estimate  of  background 
gas  velocities  falls  orders  of  magnitude  short  of  explaining  the  cooling  rate,  so  a  mechanism  is  needed 
to  convert  curl-free  Quid  expansion  which  cannot  mix  into  persistent  vorticity  which  can. 

The  symmetric  expansion  of  a  cylindrical  beam-heated  channel  causes  no  convective  mixing  per 
se.  but  asymmetries,  appearing  either  in  the  beam  deposition  or  as  density  variations  in  the  background 
gas.  permit  the  generation  of  vortex  filaments.  The  sheared  flow  attending  the  motion  and  interaction 
of  these  filaments  mixes  cold  background  gas  into  the  heated  channel  at  a  rate  which  depends  on  the 
strength  of  the  induced  vorticity  and  thus  ultimately  on  the  fluid-dynamic  asymmetries.  Three  types  of 
asymmetries  can  be  identified: 
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1.  2D  distortions  of  the  channel  from  a  circular  cross  section, 

2.  2D  asymmetries  from  pulse  displacement  off  the  channel  axis,  and 

3.  3D  distortions  such  as  characterize  a  lightning  or  spark  channel. 


In  this  paper  we  analyze  the  convective  flow  and  calculate  the  related  mixing  rate  for  the  second 
possibility,  a  lateral  displacement  of  the  second  and  subsequent  pulses  from  the  center  of  the  hot  chan¬ 
nel  created  by  the  first  pulse.  We  then  simulate  the  process  numerically  to  calibrate  our  analytic  results. 
Figure  2  provides  a  schematic  diagram  of  the  situation  which  we  are  considering.  A  hot  circular  chan¬ 
nel  of  characteristic  radius  50  and  depth  p„/po  is  assumed  to  have  reached  a  hydrostatic  equilibrium  by 
the  time  a  second  pulse  comes  along  displaced  to  the  left  by  a  distance  X&  The  initial  characteristic 
radius  of  the  second  pulse  channel  is  R0  and  it  expands  to  Ri,  eventually  coming  to  rest  there  as  a 
weak  acoustic  signal  continues  to  propagate  off  to  infinity.  This  simple  2D  situation,  shown  in  cross 
section  in  Fig.  2,  arises  when  the  second  pulse  is  being  driven  or  is  drifting  out  of  the  preformed  hot, 
low  density  channel.  Perfectly  passive  effects  such  as  transverse  flows  also  can  cause  this  displacement. 


The  radial  expansion  of  the  second  pulse  channel  generates  vorticity  according  to  the  two- 
dimensional  vorticity  equation. 


%+z-w- 

of  . 


(1) 


since  the  pressure  gradient  is  not  aligned  with  the  density  gradients  present  after  the  first  pulse.  As 
depicted  in  Fig.  3,  the  resulting  velocity  distribution  will  be  approximately  the  same  as  that  of  an 
extended  vortex  filament  pair  located  at  (3c,  ±  y)  and  having  respective  strengths  ±  K  where 


VU)  ■■  dy  J_m  dx  f .  (x,  y.  t ).  (2) 

The  bar  over  capital,  V,  has  units  of  cm2/sec  and  is  an  area  integrated  vorticity,  to  be  distinguished 

from  lower  case  v  and  v,  which  are  velocities  and  have  units  of  cm/sec.  Because  of  the  two- 
dimensional  Cartesian  geometry,  only  the  z  component  of  vorticity,  f:,  need  be  considered. 
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The  vortex  pair  shown  in  Fig.  3  propagates  horizontally  at  a  speed  which  is  roughly  1/4  that  of  the 
fluid  jetting  along  the  x  axis  between  the  two  oppositely  routing  filaments.  The  time  it  takes  a  fluid 
element  at  one  side  of  the  original  channel  (say  x  -  So)  to  reach  the  other  side  (x  -  -So)  will  be  used 
as  a  measure  of  the  mixing  time  rmix.  The  resulting  mixing  time  formula  coupled  with  our  analytic 
approximation  for  the  strength  of  the  vortex  filaments  formed,  gives  a  simple  yet  quite  accurate  expres¬ 
sion  for  scaling  laws  and  quantitative  estimates. 

In  the  next  section  the  mixing  time  formula  will  be  derived  by  analyzing  the  system  shown  in  Pig. 
3.  This  formula  is  then  coupled  with  a  dimensional-analysis  approximation  of  the  expansion-induced 
vorticity  to  show  that  mixing  of  cool  background  gas  with  the  heated  channel  interior  due  to  this  vorti- 
city  occurs  on  short  timescales,  relative  to  the  cooling  from  molecular  thermal  conduction.  The  third 
section  analyzes  the  configuration  of  Fig.  2  and  arrives  at  a  useful  analytic  expression  for  the  induced 
vorticity  to  replace  the  less  accurate  dimensional  analysis  for  K(r)  used  to  obtain  the  preliminary  esti¬ 
mates  in  Section  II.  This  analysis  gives  variations  of  integrated  vorticity  V  with  the  non-dimensional 
parameters  characterizing  the  pulse  strength,  the  lateral  pulse  displacement,  the  relative  channel  sizes, 
and  the  original  heated  channel  depth,  all  of  which  cannot  be  deduced  from  dimensional  considerations 
alone. 

Section  IV  presents  the  results  of  detailed  two-dimensional  simulations  of  the  problem  depicted 
schematically  in  Fig.  2.  These  computations  not  only  validate  the  qualitative  features  predicted  by  the 
analytic  model,  but  they  also  allow  us  to  calibrate  the  total  vorticity  integral  quantitatively  for  strong 
pulses.  This  allows  us  to  define  a  meaningful  channel  expansion  velocity  needed  for  the  theory.  Appli¬ 
cations  of  these  calibrated  formulae  to  NRL  and  LLL  data  in  Section  V  indicate  as  good  agreement 
between  predicted  mixing  times  and  observed  mixing  times  (typically  within  a  factor  of  two)  as  can-rea- 


sonably  be  expected. 
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II.  MIXING-TIME  ESTIMATE 


We  assume  that  the  asymmetry  induced  flow  field  may  be  represented  approximately  as  two  com¬ 
pact  vortex  filaments  of  strength  ±  V  at  ±y.  The  azimuthal  velocity  induced  by  each  of  these  decays  as 
1/r  away  from  the  vortex  center.  The  velocity  of  the  fluid  along  the  symmetry  plane  is  the  sum  of 
flows  induced  by  each  vortex  separately.  Figure  3  is  a  schematic  of  the  situation  with  a  sketch  of  the 
variation  of  the  fluid  flow  velocity  along  the  x-axis  shown.  This  velocity  is  given  by 


v*(x,  0)  -  2v,  cos#  -  —  -  —  — ■ — 

_  »  Vx*7F  *  X 2 +  y2 

where  cos 9  3  y/Vx2  +  y2.  The  filaments  are  migrating  in  the  same  direction  at  a  slower  velocity. 


(3) 


v,(0.J)  -  (4)- 

4 vy 

which  we  will  ignore  in  the  following  integral  estimates.  This  overall  migration  velocity  will  become 
important  when  a  series  of  pulses  is  being  considered.  In  that  case  a  quasi-steady  state  developes  in 
which  the  systematic  migration  of  the  integrated  vorticity  entrains  cold  fluid  stochastically  at  the  edges 
of  the  hot  turbulent  channel. 


To  estimate  the  mixing  time,  we  use  Eq.  (3)  to  calculate  the  time  it  takes  a  fluid  element  starting 
at  x  -  -S0  to  reach  x  —  So-  The  equation  of  motion  for  this  element  is 


i*  -  v  (x  0)  -  — _ * _ 

dt  V*U’W  w  (X2  +  y2) 


which  can  be  integrated  to  give 


So  dt  -  -  Jy  SX  +  *2>- 


'0  -  V  u  *'-•>0 

The  quadrature  can  be  performed  and  gives  the  following  mixing  time  estimate: 


(5) 


(6) 


2ir  So 

s$ 

\y\  y 

3 

(7) 


where  the  vortex  strength  V  is  given  by  Eq.  (2)  and  the  "average”  displacement  of  this  vorticity  above 
the  y  ■  0  plane  is  given  formally  by 


y  =  cfy  f-"  dx  y  e;<x,y). 


(8) 
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To  complete  this  estimate  we  need  expressions  for  V  and  y  generated  by  the  off-center  expansion 
of  the  second  pulse  channel.  The  remainder  of  this  paper  is  concerned  with  calculating  analytically  and 
calibrating  numerically  an  expression  for  V.  A  simple  dimensional  analysis  allows  us  to  make  at  least 
crude  estimates  rather  directly.  From  Eq.  (1)  it  would  appear  that  a  number  of  size  scales  enter  the 
problem:  S0,  the  radius  of  the  channel  created  by  the  first  pulse;  the  characteristic  scale  lengths  for  the 
pressure  and  density  gradients;  and  the  initial  and  final  radii  R 0  and  R  i  of  the  second,  displaced  channel 
being  formed.  Fortunately  most  of  these  scales  are  either  unimportant  or  expressible  in  terms  of  50.  If 
the  first  and  second  pulse  deposit  a  like  amount  of  energy,  R  i  and  50  will  be  roughly  equal  and  Ro  will 
be  a  modest  fraction  of  Rx.  In  the  numerical  simulations  (Fig.  7),  the  second  pulse  initial  radius 
R0  —  0.4  cm,  S0  -  1.0  cm,  and  the  overpressure  factor  for  the  first  pulse  was  30.  In  this  situation  pres¬ 
sure  balance  is  achieved  (adiabatically)  when  the  density  drops  by  a  factor  of  11.6  and  the  radius  R0 
increases  to  Rt  —  1.4  cm  —  S0. 

Other  simplifications  are  expected.  Whether  a  shock  expansion  or  an  adiabatic  expansion  is  being 
driven  by  the  energy  deposition,  the  larger  the  pressure  gradient  the  smaller  the  time  over  which  it  acts. 
Thus  the  integrated  vorticity  is  relatively  insensitive  to  the  shock  thickness.  Similarly,  the  density  gra¬ 
dient  is  integrated  over  space  so  the  inner  and  outer  densities  enter,  p0  and  p„,  but  the  scale  length  of 
the  transition  region  (shown  in  Fig.  6  as  a  band  centered  at  S  -  50)  can  be  neglected.  It  is  reasonable 
to  expect  and  will  in  fact  be  shown  that  the  maximum  vorticity  generation  occurs  when  T0  —  So-  Thus 
there  are  no  small  parameters  arising  as  ratios  of  characteristic  lengths  unless  we  consider  exceptionally 
tightly  focussed  pulses  in  which  case  R0  «  R\  —  So  or  only  slight  departures  from  the  symmetric 
superposition  of  the  pulses  in  which  case  Xo  «  So- 

The  integrated  vorticity  V  has  units  of  cm2/sec,  a  length  times  a  velocity  or  a  characteristic  time 
multiplied  by  a  velocity  squared.  The  characteristic  velocity  will  be  the  expansion  velocity.  When  the 
energy  deposition  is  fast,  as  in  the  problem  of  interest  here,  this  velocity  is  a  characteristic  sound  speed 
c,.  The  characteristic  time  is  the  expansion  duration  (if  the  energy  deposition  is  slow)  or  the  sonic 
transit  time  S<Jc,.  Thus  we  expect 
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K-c.So/*  (9) 

where  F  is  a  dimensionless  form  factor  containing  geometric  effects,  detailed  hydrodynamics  interac¬ 
tions,  and  information  about  the  channel  and  beam  profiles.  Cancellation  effects  leaving  v  •  V  v  as  the 
dominant  contribution  to  the  vorticity  source  term  generally  reduce  F  somewhat  below  unity.  This  will 
be  demonstrated  below. 


If  we  assume  that  y 


S0,  Eq.  (7)  for  the  mixing  time  becomes 


7.7— . 

For  a  strong  pulse  which  heats  the  gas  appreciably  on  passage,  we  shall  assume  F 


(10) 


— .  Taking 


5  x  104  cm/sec  as  a  generic  value  for  c,  and  choosing  S0  —  1  cm  gives 

7.7  x  1  __ 

Tmix  ~  I  ■  300  n  sec. 

5x  lf^xi 

This  300  /usee  'mixing  time”  from  the  dimensional  analysis  estimate  of  Eq.  (10)  is  the  time  for  cold 
material  from  one  side  of  a  hot  channel  to  cross  the  channel  and  reach  the  other  side.  True  mixing  of 
hot  and  cold  fluid  will  probably  take  somewhat  longer  because  one  or  two  rotations  of  the  vortices  will 
be  required  to  entrain  and  smear  in  an  appreciable  amount  of  the  colder  fluid.  This  'mixing'  time, 
however  it  is  best  defined,  is  clearly  important  when  comparable  to  only  a  few  pulse  separation  times. 


In  the  next  section  an  integration  of  the  vorticity  source  term  in  Eq.  (1)  is  performed  to  improve 
on  Eqs.  (9)  and  (10),  our  dimensional  analysis  estimates.  Using  a  few  reasonable  assumptions,  a  qua¬ 
drature  for  integrated  vorticity  V  is  obtained  which  will  display  the  various  non-dimensional  depen¬ 
dences  of  the  form/fudge  factor  Fin  Eqs.  (9)  and  (10). 


III.  AN  APPROXIMATE  VORTICITY  INTEGRAL 

Equation  (2)  defines  the  integrated  vortex  strength  at  any  instant  in  time  and  will  be  used  in  con¬ 
junction  with  Eq.  (1)  to  develop  an  approximate  analytic  expression  to  replace  Eq.  (9)  in  the  mixing 
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time  estimate  Eq.  (7).  For  our  purposes  we  only  need  V(t)%  where  r  is  a  time  sufficiently  long  that 
expansion  and  vorticity  generation  is  complete.  Using  the  geometry  of  Fig.  1, 


~  So  dt  So' d9  f*u>rdr 


Vp  x  V/> 


-H, 


(id 


where  the  spatial  integral  is  performed  over  the  fluid  outside  the  expanding  cylinder  R0  <  R  (/)  <  Rx. 


We  assume  that  this  outside  region  ( r  >  R(i ))  behaves  incompressibly  and  the  inside,  beam- 
heated  region  expands  uniformly.  Thus 

v,(r,r)  -  for  r  <  (r)  and  (12) 

v,(r,/)  -  - for  r  >  R(t) 

r 

specifies  the  flow  everywhere  as  a  function  of  the  heated  region  radius  R  (r)  and  velocity  U(t)  =  R  (t). 
This  flow  has  a  uniform  but  time-varying  divergence  inside  the  heated  region  and  zero  divergence  out¬ 
side.  In  fact,  the  fluid  inside  the  just-heated  channel  will  give  up  energy  to  the  cold  surrounding  fluid 
via  shocks  and  a  fraction  of  the  pulse-deposited  energy  will  even  escape  to  °°  as  an  acoustic  wave.  The 
smooth  shape  of  the  expanded  channels  observed  from  ID  hole-boring  calculations3  and  the  computa¬ 
tional  fact  that  most  of  the  deposited  energy  stays  close  to  the  original  pulse  deposition  lends  hope  for 
the  approximations  implied  by  using  the  flow  field  (12). 


As  will  soon  be  shown  somewhat  more  explicitly,  the  important  feature  of  vorticity  generation  is 
how  far  radially  each  fluid  element  moves.  Pressure  gradients  arising  from  accelerations  of  a  fluid  ele¬ 
ment  do  not  really  contribute  to  vorticity  in  the  present  context.  Because  the  fluid  elements  begin  and 
end  their  expansion-induced  displacement  at  rest,  the  average  acceleration  is  zero.  The  v-  V  v  vorticity 
source  term  has  the  same  sign  throughout  the  expansion  and  consequently  contributes  more  strongly  to 


the  integrated  vorticity. 
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The  velocity  field  (12)  can  be  substituted  into  the  fluid  momentum  equation  for  the  outside 
region  and  yields 


d\r  9v,  |  fly.  1  d(UR)  U2R 1 

dt  dt  V'  dr  r  dt  r3 


Since  we  are  going  to  assume  the  motion  behaves  incompressibly  in  the  outer  region  r  >  R{t ),  the 
second  term  on  the  right  hand  side  of  Eq.  (11)  vanishes  simplifying  the  integral  considerably.  The 
acceleration  is  being  taken  in  the  radial  direction  with  values  strictly  appropriate  only  when  feedback  of 
the  asymmetric  density  gradients  on  the  driving  expansion  flow  are  small. 


The  theoretical  model  appears  to  work  quite  well  for  large  density  variations  as  well,  a  result  that 
is  understandable  in  hindsight.  The  maximum  vortex-induced  flow  speed  which  resides  in  the  system  is 
much  smaller  than  the  maximum  expansion  speed.  The  vorticity  is  generated  essentially  instantane¬ 
ously  relative  to  the  mixing  timescale  r^.  Thus  the  generation  term  can  be  calculated  assuming  that 
the  initial  asymmetric  density  gradients  do  not  change  the  expansion-driven  pressure  gradients  and  that 
the  vorticity  which  develops  does  not  affect  the  density  gradients  during  the  relatively  brief  expansion. 


Under  these  assumptions  only  the  density  gradient  in  the  9  direction  need  be  crossed  with  the 
radial  pressure  gradients.  Using  Eq.  (13)  and 


in  Eq.  (11)  we  obtain 


1  dp 

7  d9 


_ 1  Bp 

(UR)1 

1  9 (UR) 

pr  90 

r3 

r  dt 

We  now  change  to  the  Lagrangian  variables  (r0,  0O)  where 


r  —  r(/)  -  yfr^~7~RHt)~-Rj  and 
0-0o- 


(14) 


(15) 


(16) 


8 


. 

,-, _ _ _  _ 
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Then 


rdr  —  r0dr0  (17) 

and  r(r)  is  the  instantaneous  radial  position  of  a  fluid  element  initially  at  r(0)  —  r0.  Because  the  flow 


described  by  Eq.  (12)  is  incompressible,  the  continuity  equation  becomes 

pV-K-0, 


ifi. _ 

dt 


and  therefore 


p (r(/),  9)  -  p(r0,  9). 


Thus  Eq.  (15)  becomes 

P(r)  -  />/>/,>, 


(UR)2 


1  9(UR) 


r3(t ) 


r(t)  dt 


(18) 


(19) 


(20) 


Note  that  the  9  gradient  of  p  decreases  as  the  expansion  proceeds  because  the  fluid  elements  are  mov¬ 
ing  to  larger  radius  while  retaining  the  same  9  variation. 

At  this  point  we  assume  a  general  form  for  the  density  profile  pis)  in  the  original  low  density 
channel. 


where 


Then 


s  —  Xtf  +  rtf  -  2X0r0  cos 9, 
g( 0)  -  1  and  g(<»)  —  0. 


(21) 


(22) 


v»p  _ 

1  dp  . 

P- 

1  dg(s) 

p 

rp  99 

P  0 

r  99 

(23) 


and  the  9  integral  in  Eq.  (15)  can  be  performed  without  approximation. 
This  gives  us 

So  d9*ht  “  g(s{r°’  ~  0)) 


giro  +  Xq)  —  giX o  —  rg)  ro  ^  Xo 
giro  +  Xq)  —  giro  —  Xq)  ro  >  Xo 


(24) 
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and  Eq.  (20)  becomes 


o  r2(t)  r2(t) 


g(rQ  +  X0) 


—  dr0r0  (UR)i 
x*  r2U)  r2(t) 


*<'o-*o>)g5) 


Still  working  on  this  rather  general  form  of  the  equation,  approximations  must  be  made  to  get  a 
useable  analytic  result  because  Eq.  (25)  is  a  hairy  2D  quadrature  even  assuming  that  the  appropriate 
R  ( t )  for  our  particular  case  is  known.  The  approximation  most  useful  is  to  evaluate  the  integrand  at 
one  time.  If  we  assume  U (r)  —  U(O)  **  0  so  that  the  initial  and  final  states  of  the  expansion  are  at 

a  — 

rest,  the  term  —  (UR)  integrates  to  zero.  To  approximate  V(r),  let  us  evaluate  /•(/)  in  the  integrand 
ot 

a 

at  the  time  t„  when  the  expansion  flux  RmUm  =  R Um)  U(tm)  is  maximum.  Then  the  —  (UR)  term 
vanishes  and  r(t)  is  still  a  function  of  r0,  the  initial  fluid  element  radius,  but  not  a  function  of  t. 


rm  *  y/r$  +  Rj,  -  R$ . 


The  desired  integral  becomes 


-  I Jr  «  +  *o) 

|poJJo  V  o  r” 


In  Eq.  (27)  the  time  integral  can  be  simply  approximated  as  U2R2  t/2  where  r  -  2(R\  -  Ro)/Um . 


ti y  ,  U2r  .  lp-\lrm.  R£g(ro  +  Xo) 

rx o  RmgiX o  ~  rp)  fm  j 

~  dr°r0  lrl  +  *2  _  D  l\2  ~  Jx0  dr°r° 


R2g(r0-X0 ) 
W  +  Rj,  -  /?o  >2 
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Now  the  integral  is  written  in  terms  of  three  non-dimensional  parameters 
a  as  Rq/ Rm  <1, 

b  =  RJS0  <  1,  and  (29) 

c  a  X<JS0  -  1. 

where  S0  is  a  characteristic  scale  length  for  variation  of  g(s/S0)  and  hence  of  the  original  channel  den¬ 
sity  profile,  Eq.  (21).  Letting  the  integration  variable  be  t>  =  r/Rm  yields 


\s: 


di\t) 


£(S0fob  +  c]) 


(1  +  n2  —  a2)2 
,  ^(Solc-Tib])  f00  ^(Sofob-  cl)  1 

'J.  <r+”a5)3  (TTV^Pi7) 


(30) 


If  the  remaining  composite  integral  is  identified  as  a  three  parameter  form  factor  /2(a,b,c)  with 
values  of  order  unity  and  if  (Jm  —  c,,  an  assignment  which  will  be  justified  empirically  by  detailed 
numerical  simulation  in  the  next  section  for  strong  pulses,  the  very  useful  improvement  on  the  dimen¬ 
sional  approximation  Eq.  (9)  is  obtained: 


K(r)  =  c,(R\ ~  *oMn 
Here  the  form  factor  /2,  defined  as 


P  o 


/2(a,b,c). 


(31) 


,  _  r"  .  g(50fc»b+cl)  fdb  g(Solc-Tjbl) 

/Afcd-J.  +  *>n  (1+v_aV 

g(S0[T)b  - 
(1  +  t}2  —  a2)2 ' 

is  generally  less  than  1/2.  The  sign  of  /2  indicates  the  direction  of  flow  in  the  .v.^plane— 


r°°  ,  g(S0[T)b- 

Jc/b^^d+^-a2)2’ 


(32) 


counterclockwise  or  clockwise. 


The  general  integral  form  factor,  Eq.  (32),  can  be  integrated  numerically  for  any  reasonable 
profile  g(s/S0).  Figure  4  shows  plots  of  /2  versus  c,  the  measure  of  channel  separation  for  several 
values  of  a  at  b  -  0.7  and  for  several  values  of  b  at  a  -  0.6,  values  closely  approximating  the  actual 
initial  conditions  for  the  detailed  calculation  of  the  next  section.  In  this  figure  the  super  Gaussian 
profile  was  used. 
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gsais/So)  -  e~,1/si  (33) 

and  the  integrations  performed  using  a  numerical  quadrature  algorithm.  Calculations  have  also  been 

carried  out  using  a  Bennett  profile. 


/sw(a,b,c)  -  J  (1  +  c)2  +  b^(l  -  a2)  (1  -  c)2  +  b*(l  -  a2)  °6) 

for  ab  <  c  -  1.  This  is  the  case  when  the  second  pulse  is  wholly  outside  the  channel  formed  by  the 
first  pulse.  As  c  approaches  fD  decreases  as  -l/c3  so  pulse  channels  separated  by  more  than  3  radii 
interact  only  weakly. 


When  the  channels  are  closer,  ab  >  c  —  1  and  the  channels  overlap  somewhat  initially.  When 
-ab  <  c  -  1  there  is  still  part  of  the  second  channel  of  initial  radius  R0  outside  the  cylinder  s  -  S0  so 


/sw(a,b,c)  —  J  (1+c)2+*b2(1_a2)  1  • 

Equation  (37)  agrees  with  (36)  at  the  interfacial  separation  c  -  1  +  ab. 


There  is  a  third  region  when  the  initial  radius  of  the  second  channel  lies  wholly  within  the  first, 
i.e.,  S0  >  X0  +  R0  or  1  >  c  +  ab.  Then  /Sw  (a,b,c)  satisfies  Eq.  (36)  again.  The  cross  over  point  at  c 
+  ab  -  1  again  has  the  two  expressions  equal.  In  the  intermediate  region  (1  -  ab  <  c  <  1  +  ab), 
/sw  takes  on  its  peak  (negative)  value  at  the  exterior  touching  point  c  ■  1  +  ab.  For  larger  values  of 
c  the  magnitude  l/$wl  decreases  and  for  smaller  values  |/Swl  decreases  monotonically  to  zero  at  c  —  0. 
This  peak  negative  value  is 
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When  ab  -  b2  -  1/2,  |/SWIB 
respectively. 


f  »  2(1  +  ab) 

/sWmM  (4  +  4ab  +  b2)  ( 

:  -  6/13  which  is  close  to  the  values  in  Figs.  4  and  5  for  /so  and  f8 


Also  of  interest  is  the  slope  at  c  -  0  where  the  pulses  are  concentric  and  no  vorticity  genera¬ 


tion  is  expected: 


d _ 

dcL  1  + 


-2b2 

tfU  -  a2)]2 


which  is  non-zero  and  large.  Thus  even  modest  non-concentricity  leads  to  appreciable  vorticity  and 
mixing.  Similar  behavior  is  observed  in  the  two  smooth  profiles  considered  as  well.  Figure  6  compares 

h  |^j-.  c  ,  /sg  c|  and  /sw  j-4j.  c|  to  demonstrate  the  similarity  among  the 

vorticity  generation  form  factors  regardless  of  the  profile.  This  relative  profile  insensitivity  arises 
because  V  is  an  integral  quantity,  the  total  upper  half  plane  circulation. 

Returning  to  our  integrated  vorticity  estimate  of  Eq.  (31)  with  a  maximum  form  factor  of  6/13, 


1^1™*=  4- c,  (*,-*o)  In  N-  • 
13  \  p  o 


Using  vieI  as  the  maximum  on  axis  velocity  and  assuming  y  ~  R{  -  R0  gives 


1Vjrlm“=S  13  it  ln(p0  la¬ 


under  optimum  conditions  the  velocity  between  the  vortex  filaments  approaches  a  quarter  of  the  sound 
speed  in  the  surrounding  fluid.  When  the  expansion  is  subsonic  throughout  (where  this  analysis  ought 
to  be  most  accurate),  c,  in  Eq.  (41)  should  be  replaced  by  the  appropriately  averaged  expansion  velocity 


Using  Eq.  (7)  for  the  mixing  time  with  y  =  S0  *=  /?i  -  /?0  and  I  l/!mtx  from  Eq.  (40) 


above  gives 
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With  c,  as  5  x  104  cm/sec,  In  (p«/po)  ~  2.5,  and  S0  —  I  cm,  the  fastest  mixing  time  is  seen  to  be 
about  160  psec.  This  estimated  fastest  mixing  time  is  within  a  factor  of  two  of  the  dimensional  esti¬ 
mate.  In  any  particular  system  the  integral  of  Eq.  (32)  can  be  performed  numerically  and  the  resulting 
form  factor  substituted  into  Eq.  (25)  for  the  integrated  vorticity. 

The  sign  of  ^(r)  is  negative  for  the  configuration  of  Fig.  2  where  the  second  pulse  is  centered  to 
the  left  of  the  original  channel.  This  means  that  the  jet  of  colder  wall  gas  across  the  original  channel 
starts  on  the  side  opposite  to  the  displaced  second  pulse  and  rushes  toward  the  newly  heated  region. 
On  average  the  old  channel  moves  toward  the  new  channel.  This  average  motion  does  not  imply  mix¬ 
ing  but  the  spatial  behavior  of  the  flow  from  the  vortex  filament  pair,  as  will  be  seen  in  the  simulations 
of  Section  IV,  effectively  bisects  the  composite  hot  channel  with  a  cold  jet.  The  result  is  two  smaller 
channels  above  and  below  the  original  symmetry  plane  in  which  the  vortex  filaments  are  now  centered. 
A  third  pulse  located  at  y  -  0  between  the  two  modified  channels  will  cause  each  of  these  smaller  chan¬ 
nels  to  bifurcate  again  with  fluid  jets  from  the  top  of  the  upper  channel  and  the  bottom  of  the  lower 
channel  impinging  on  the  third  expanding  channel  from  above  and  below. 

In  the  next  section  the  problem  of  Fig.  7  is  solved  numerically  with  p„/p 0  -  10,  R0  -  0.4  cm, 
R  |  —  1.4  cm,  S0  —  1.0  cm,  and  X0  —  1.0  cm.  If  we  assume  R,„  —  1.0  cm,  a  —  0.4,  b  —  1.0,  c  •  1.0, 
then  /j(.4. 1.0, 1.0)  =  -1/5  from  Fig.  16.  Substituting  these  number*  into  Eq.  (31)  with  c,  taken  as 
5  x  104  cm/sec  gives  V(r)  =  -2.3  xlO4  cmJ/sec.  This  number  is  within  10%  of  the  measured  value  - 
2.1  x  104  cm2/sec  (Section  IV).  The  value  of  c,  used  should  be  evaluated  in  the  original  channel  den¬ 
sity  and  temperature  profiles  at  the  location  of  the  center  of  the  second  pulse,  i.e.,  x  -  -  X0  -  -1  cm 
for  the  example  here.  The  characteristic  time  r  must  be  equal  to  (/?  t  -  rt0)/(c,/2)  —25  psec  to  justify 
using  c,  as  the  effective  maximum  expansion  velocity  of  the  heated  cylinder.  This  value  is  also  within 
10%  of  the  observed  value  of  r  <  30  psec. 
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Using  Eq.  (3)  with  x  -  0  gives  the  maximum  flow  velocity  expected  in  the  symmetry  plane, 

v,  -  -  =  -90  m/sec,  (43) 

-try 

when  y  —  V!  cm.  This  value  is  in  good  agreement  with  the  measured  jet  velocity  of  -85  m/sec 

seen  in  the  detailed  simulations.  The  mixing  time  calculated  from  Eq.  (7)  for  these  parameters  is  about 
330  fisec. 

IV.  DETAILED  NUMERICAL  SIMULATIONS 

Preliminary  detailed  numerical  simulations  of  this  problem  have  been  performed  using  the 
FAST2D4-5  6  computer  code  to  validate  and  calibrate  the  approximate  analytic  model  developed  above. 
In  addition  to  accounting  for  shocks  properly,  which  the  theoretical  model  does  not  do,  the  simulations 
are  capable  of  describing  the  late-time  motions  and  profiles  as  modified  by  the  induced  vorticity.  Thus 
the  various  elements  of  the  mixing- time  approximations,  and  their  significance,  can  be  evaluated  using 
detailed  calculations  of  representative  problems.  In  this  paper  we  restrict  ourselves  to  a  few  calculations 
related  to  the  specific  configuration  of  Fig.  7  to  calibrate  and  validate  the  most  questionable  assump¬ 
tions  leading  to  Eq.  (28).  Questions  of  channel  cooling  and  complex  multipulse  flows  are  left  for  future 
studies. 

The  calculations  based  on  fig.  7  used  Cartesian  coordinates  on  a  stationary  rectangular  grid  con¬ 
sisting  of  100  by  30  cells  of  dimensions  8x  and  By,  each  varying  from  1  mm  to  5  mm.  As  shown  in 
Figs.  8  and  9,  an  embedded,  finely  zoned  region  of  50  x  20  cells  1  mm  on  a  side  was  centered  on  the 
initial  1.0  cm  radius  hot  channel.  With  increasing  distance  from  the  center  of  the  channel,  the  values 
of  8x  and  8 y  transitioned  smoothly  from  1  mm  to  5  mm.  Near  the  sides  and  top  of  the  grid,  Sx  and  By 
were  both  5  mm. 

Figure  10  shows  a  plot  of  p  (x,  0, 0)  and  Pix,  0,0)  along  the  y  —  0  plane  at  the  time  the  second 
pulse  is  deposited.  Bennett  profiles  l/tl  +  s2/S$  ]2  for  the  density  and  pressure  deviation  from  ambient 
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were  used  in  these  calculations.  The  dots  (•)  on  the  two  curves  show  the  location  and  spacing  (1  mm) 
of  the  finite  difference  grid  points  in  the  vicinity  of  the  channel  centers. 

To  compute  the  vorticity  approximately  we  used  the  relation 
Y(t)  -  SI*  t  dx  dy  -  \x(x,  O.r)  dx 

a  f  "**  vx(x,0.t)dx  (44) 

where  and  are  the  respective  limits  of  the  coordinates  on  our  grid.  To  derive  Eq.  (44)  we 
have  assumed  that  vx,  vy  —  0  as  x,y  —  <». 

Figure  11  shows  the  results  of  our  calculations  of  Y(t)  for  R0  -  0.4,  S0  -  1.0,  and  X0  —  1.0  cm 
using  the  density  profile  of  the  channel  formed  by  the  first  pulse  and  the  pressure  profile  of  the  second 
pulse  as  shown  in  Fig.  10.  After  an  initial  positive  transient,  the  integrated  vorticity  settles  down  to  a 
steady  state  negative  value  after  a  weak  relaxation  oscillation. 

dv 

In  deriving  Eq.  (27)  from  Eq.  (25)  the  time  derivative  term  -r~,  which  must  integrate  to  zero 

ot 

dv 

during  the  expansion,  is  neglected  in  favor  of  the  v  —  term  which  is  always  of  the  same  sign  and 

or 

hence  contributes  most  to  the  integral.  A  posteriori  justification  for  this  approximation  is  obtained 
from  Fig.  11  which  shows  the  general  cancellation  of  a  big  transient  component  in  Yft),  as  assumed  in 
Section  III,  and  good  quantitative  agreement  with  the  residual  secular  term.  This  means  that  the 
amount  of  vorticity  generated  depends  strongly  on  the  eventual  radial  displacement  of  each  fluid  ele¬ 
ment  but  not  on  the  detailed  expansion  history.  The  scale  of  vorticity  is  set  by  an  average  expansion 
velocity  only. 


Figure  12  shows  a  schematic  diagram  to  explain  physically  how  the  initial  positive  vorticity  spike 
of  Fig.  1 1  occurs.  The  expanding  shock,  shown  at  four  different  times  as  a  dashed  line,  travels  faster 
and  drives  the  post  shock  flow  faster  in  the  low  density  material.  The  density  gradient  is  shown  here 
localized  between  the  two  solid  semi-circles  in  an  annulus  centered  at  50.  The  shear  resulting  from  the 
differential  acceleration  of  the  inner  and  the  outer  fluid  corresponds  to  vorticity  of  positive  sign  which 
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vanishes,  as  assumed,  when  differential  velocity  ceases.  The  angle  between  Vp  and  V/*  (assumed  nor¬ 
mal  to  the  dashed  line  shock  fronts)  is  shown  at  three  locations  along  the  original  channel  periphery. 
The  residual  vorticity,  which  corresponds  to  that  appearing  in  Eqs.  (27)  and  (30),  in  effect  measures 
the  net  outward  expansion  of  the  new  channel. 

Figures  13  and  14  show  the  vorticity  strength  F(r)  as  a  function  of  X0  with  the  other  remaining 
variables  in  our  standard  case  held  constant.  When  the  axis  of  the  second  pulse  fell  within  50  of  the 
center  of  the  channel  (Fig.  13),  the  transient  term  appeared  to  oscillate  more  rapidly  than  the  cases  for 
which  the  axis  of  the  second  pulse  fell  ’outside”  the  original  channel  (Fig.  14).  The  residual  values 
V(t)  to  which  the  integrated  vorticities  in  Fig.  13  related  were  also  more  nearly  equal  than  in  Fig.  14. 
The  general  trend  of  VU  —  <»)  in  the  five  cases  of  Figs.  13  and  14  follows  qualitatively  the  shape  of 
/2(a,b,c)  plotted  for  various  cases  in  Figs.  4,  S,  and  6.  The  uncancelled  vorticity,  which  we  called  F(r) 
in  Section  III,  decreases  from  a  maximum  around  X0  —  S 0  both  as  X0  — ’  0  and  as  X0  gets  large  com¬ 
pared  to  S0-  The  ratio  of  V(r)  at  Xq  -  0.2  cm  to  K(r)  at  X0  ■»  0.1  cm  is  probably  not  as  small  as 
expected  from  the  curves  of  /2  in  Fig.  5  because  the  sound  speed  is  higher  when  the  second  pulse  is 
near  the  center  of  the  first  channel.  This  increasing  sound  speed  partially  counterbalances  the  tendency 
of  /2  to  fall  off  as  X0  approaches  zero. 

The  major  issue  in  interpreting  these  calculations  using  the  formulae  of  Section  III  revolves 

Uh 

around  the  choice  of  — --  for  Eq.  (30).  For  Eq.  (31)  and  later  estimates  we  have  assumed  Um  is  the 

sound  speed  of  the  fluid  at  the  second  pulse  center  before  the  second  pulse  is  laid  down.  We  have  also 
assumed  that  the  average  expansion  velocity  is  UJ  2. 

Figure  IS  compares  the  effects  of  initial  pulse  size  on  the  second  pulse  with  the  same  values  of 
the  channel  parameters  and  the  overpressure  and  displacement  X0  of  the  second  pulse.  We  notice  that 
the  transient  term  appears  to  scale  in  magnitude  as  R0  while  the  value  of  the  "residual”  vorticity  appears 
to  scale  as  R$.  We  may  explain  the  scaling  of  the  residual  vorticity  in  terms  of  our  analytic  "late  time" 
treatment.  Instead  of  using  the  estimate 
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Um  -  C„ 

as  in  Eq.  (31),  we  may  calculate  Um  from  the  equation 


(45) 


Um 


2(/t,  —  /Ip) 

T 


(46) 


where  r  is  the  time  of  expansion  from  A0  to  R\.  In  Fig.  15,  we  see  that  the  extrema  in  V  for  the  cases 
R0  “  0.2  cm  and  0.4  cm  occur  at  almost  the  same  values  of  t  and  that  the  relaxation  to  the  final  value 
of  V  occurs  at  the  same  rate  for  the  two  cases.  Thus  we  have  (for  any  consistent  definition  of  r)  the 
relation 


However  we  expect  that 


a  t*. 


since 


Thus  Eq.  (46)  yields 


Estimates  of  the  form  factor  indicate  that 


R\a 


\b 


V*a  ~  \umB. 


fu  *  flB- 

Using  Eqs.  (31)  and  (46H51),  we  find  that 

W  ,  ( Ru-Rm )2  ^  1 
VB(T)  ( R\b  ~  Rob)*  ^ 

for  the  case  shown  in  Fig.  15. 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 


This  analysis  indicates  that  Eq.  (45)  provides  only  a  crude  estimate  of  Um  as  the  energy  deposited 
becomes  small.  Oearly  the  relaxation  time  of  the  transients  depends  as  much  on  the  already  existing 
density  gradients  as  it  does  on  the  new  channel  expansion  laid  down.  The  ambient  sound  speed  at 
X0  -  1.0  cm  is  about  4  x  104  cm/sec  for  the  standard  case  and  the  reduced  energy  case  of  Fig.  15.  The 
final  relaxation  to  the  long-time  value  V(r)  begins  at  about  40  usee  for  both  cases,  about  the  time  it 
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takes  a  sound  wave  to  reach  the  far  side  of  the  original  hot  channel  from  the  edge  of  the  second  pulse 
(a  distance  of  1.8  cm  and  1.6  cm  for  cases  A  and  B  respectively). 

Clearly  additional  work  is  required  to  refine  the  approximation  of  Eq.  (31),  to  improve  the  rather 
cavalier  way  in  which  Um  has  been  determined.  We  expect  to  carry  this  study  forward  in  a  future 
report  using  cases  for  analysis  and  simulation  which  are  of  prime  experimental  interest. 

In  Table  I  and  Fig.  16,  we  present  estimates  of  the  form  factor  /2  for  the  case  So  ~  1.0  cm, 
R0  -  0.4  cm,  R\  -  1.4  cm,  p«/po  -  10  (channel  formed  by  first  pulse)  and  PJPQ  -  31  (second 
pulse),  as  shown  in  Figs.  7  and  10.  These  values  correspond  to  a=  0.4  and  ba  1.0.  We  vary  c  from 
0.2  to  2.0.  To  estimate  the  value  of  £/m,  we  used  Eqs.  (45)  and  (46).  For  Eq.  (46),  we  define  r  to  be 
the  time  elapsed  between  the  deposition  of  the  second  pulse  (r  -  0  in  Figs.  13  and  14)  and  the  last 
extremum  in  V(t)  before  the  monotonic  relaxation  to  the  residual  value.  The  solid  curve  is  the  esti¬ 
mate  from  Eq.  (32). 

For  both  methods  of  estimating  U„  the  values  of  fi  appear  to  be  somewhat  low  near  c  —  1.0,  as 
we  see  from  comparison  of  the  simulation  data  points  with  the  analytic  calculations  also  shown  in  Fig. 
16  for  the  Bennett  profile,  which  was  used  in  our  simulations.  Here  Eq.  (45)  appears  to  provide  a 
better  estimate  of  Um  than  does  Eq.  (46),  at  least  when  R0  -  1.0  cm.  Better  agreement  is  expected 
when  the  energy  deposition  is  slow  and  the  adiabatic  expansion  treatment  is  a  better  approximation. 

Finally,  returning  to  Fig.  1,  assume  Um  =  2.0  x  104  cm/ sec,  R\  -  R0  =  15%  of  4  mm  ”  0.06 
cm.  In  (p«/po)  ~  0.33,  and  choose  f\  -  -1/20  for  a  smaller  relative  displacement  than  shown  in  Fig. 
3.  The  expansion  velocity  is  smaller  here  because  the  pulse  is  weak.  The  predicted  value  of  K(t)  is  20 
cm2/ sec.  The  corresponding  mixing  time  is  about  60  msec  for  two  pulses  and  perhaps  a  half  or  a  third 
of  that  time  if  the  vorticity  from  the  subsequent  third  and  fourth  pulses  add  roughly  linearly.  More 
accurate  estimates  for  these  data,  kindly  provided  by  the  authors,2  must  await  specific  detailed  calcula¬ 
tions  in  a  future  publication. 
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V.  SUMMARY  AND  CONCLUSIONS 


Formulae  have  been  developed  and  calibrated  for  calculating  the  integrated  vorticity  and  subse¬ 
quent  convective  mixing  time  resulting  after  a  second  circular  heating  pulse  is  injected  parallel  to  the 
axis  but  off  center  in  a  hot  channel  created  by  an  earlier  pulse.  The  vorticity  arises  because  the  density 
gradient  asymmetries  interact  with  pressure  gradients  attendant  on  the  second  heating  pulse.  Two 
oppositely  routing  vortex  filaments  give  rise  to  a  shearing,  mixing  flow  field  which  tends  to  convect  the 
original  channel  toward  the  center  of  the  second  pulse. 


The  following  formulae  summarize  the  results  of  this  paper.  Assume  a  hot  channel  exists  with  a 
central  density  p0  lower  than  ambient  density  pM  and  with  a  characteristic  radius  S0  for  the  density  vari¬ 
ation.  In  the  simulations  performed  for  this  paper  a  Bennett  profile  was  used: 


P(S)  “  Poo  ~ 


(p<»  -  Po) 


(53) 


11  +  (s/So)2!2 

both  for  the  original  channel  and  for  the  subsequent  energy  depositions.  A  second  pulse  is  deposited  a 
distance  X0  off  center  in  the  heated  channel  with  an  overpressure  8 /*„,,*  at  the  center  and  a  characteris¬ 
tic  radius  Ro- 


P(r) 


Poo  + 


8  P„ 


[I  +  (r/R0)2]2 


(54) 


for  the  simulations  reported  here.  The  extra  energy  deposited  per  unit  length  by  this  second  pulse  is 


IE 


ItRp  8P mm 

(y  —  1) 


which  may  also  be  used  to  calculate  from  A£.  For  the  calculations  reported  here  8  Pm,, 


(55) 
30  and 


P0  *  0.4  cm  (generally). 


If  adiabatic  expansion  of  the  second  pulse  channel  to  ambient  pressure  is  an  adequate  approxima¬ 
tion  as  suggested  by  Lampe,3  7  the  final  characteristic  radius  of  the  expanded  pulse  channel  will  be 
(approximately) 
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*  « 


Giving  R\  —  Rq  for  use  in  Eq.  (31).  The  integrated  vorticity  in  each  of  the  two  filaments,  ±  V(t),  is 
given  adequately  by  the  formula 


p(r)  *  cs(X0HRi  -  Ro)  in  —  h  4L’  (57) 

l  PO  Rm  ^0  ^0 

where  the  sound  speed  c,  (XQ)  is  evaluated  in  the  original  channel  at  the  location  where  the  second 
pulse  will  be  centered.  The  radius  Rm  lies  between  Ro  and  R\  and  we  have  taken  it  to  be  1.0  cm  when 
R0m.  o.4  cm  and  -  30.  If  the  original  channel  overpressure  is  S/^,  p„/po  is  given  by  an 
expression  analogous  to  Eq.  (36) 


wpo>- 

*  oo 


+  i  r. 


The  geometric  form  factor  f2  (a,b,c)  is  given  by  the  simple  quadrature 


h  (a,b,c)  -  r\di\ 


tglSpfab  -1-  c])  —  glSofob  —  cl)] 
[I  +  t)2  -  a2]2 


where  the  original  channel  density  satisfies 


-In  —  *<j/S0> 

p(s)-p  »e  lPo1  .  (60) 

Equation  (59)  for  the  form  factor  is  a  permissible  simplification  when  g{s/S0)  is  expressible  in  even 

powers  of  s.  Plots  of  this  form  factor  for  a  number  of  cases  are  presented  in  Figs.  4,  5,  6  and  16.  It 
typically  takes  on  negative  values  with  an  extremum  of  f2  —  y  in  the  vicinity  of  X0  *  So  (c  —  1)  and 
approaches  zero  for  both  large  and  small  values  of  X0. 

Using  Eq.  (57)  for  the  generated  vorticity,  a  mixing  time  estimate  can  be  derived  as  in  Section  II, 


>/6ir  Sj 


Tmi*~  |K(r)| 
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This  is  the  time  it  takes  coid  fluid  from  one  side  of  the  original  channel  to  cross  the  hot  region  and 
reach  the  other  side— thus  effectively  bisecting  the  hot  channel  with  cool  fluid.  The  negative  sign  on  /2 
means  physically  that  the  vortex-induced  cold  jet  moves  toward  and  into  the  second  pulse  channel 
rather  than  away  from  it. 


As  an  example  of  the  use  of  these  formulae,  consider  the  case  of  A£  ergs/cm  being  deposited  by 
each  pulse  in  a  cylinder  of  initial  radius  Ro ■  The  overpressure  factor  is 

_ _  Syntax  (y  —  1)  A£ 

r  =  —  -  m 
and  the  mixing  time  estimate  of  Eq.  (61)  can  then  be  written 


i 


r  o _ (l  +  ny _ 

TmiX  Ur  -L 

1(1  +  D 2r  -  l]ln(l  +  r> 

The  typical  turbulence  velocity  UT  is  given  by  the  expression 


(63) 


UT  -  ®  103  cm/sec  (64) 

where  the  value  of  103  cm/sec  assumes  <  /2>c,  —  104cm/sec.  Since  cs  —  5  x  104  cm/sec  and 

|/2|  —  j  in  the  range  0.5  S0  <  X0  <  1.5  S0,  this  is  a  good  approximation. 


A  single  curve,  as  seen  from  Eq.  (63),  serves  to  characterize  the  non-dimensional  vorticity  mixing 
time  in  terms  of  the  instantaneous  pulse  overpressure.  However,  since  8PmjP*>  also  depends  on  /?0, 
curves  of  rmjx  for  several  values  of  Ro  are  plotted  in  Fig.  17.  For  tightly  focused  beams  (£0  small  and 
£  large)  the  mixing  times  are  very  fast  so  convective/turbulent  cooling  is  bound  to  play  a  significant 
rote  in  muitipulse  propagation.  Overpressures  significantly  less  than  unity  must  be  used  with  pulses 
separated  by  10~3  sec  or  more  when  the  beam  radius  is  1  mm.  In  fact  mixing  will  significantly  alter 
channel  dynamics  with  a  deposition  radius  of  1  mm  and  overpressure  of  10  unless  the  pulse  separation 
in  time  is  reduced  well  below  100  ^isec. 
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There  is  much  still  to  be  learned  about  this  vorticity  generation  mechanism  which  may  have  a  tel¬ 
ling  role  in  determining  overall  channel  coding  characteristics.  In  closing  this  paper  it  is  appropriate  to 
identify  several  of  these  issues  which  appear  to  play  a  major  role  and  which  have  been  dealt  with  only 
approximately  above.  Where  possible  the  likely  effect  is  described  and  discussed. 

In  our  theoretical  model  for  the  integrated  vorticity.  a  great  simplification  was  achieved  by 
evaluating  the  integrand  at  a  point  in  time  when  the  maximum  expansion  velocity  was  attained  and 

then  replacing  Ul  y  with  cs(Rt  -  R0).  For  moderate  strength  shocks  this  should  work  reasonably 

well.  Ignoring  for  now  the  question  of  what  is  the  "correct"  average  to  take,  it  should  be  clear  that  the 
formula  is  probably  wrong  for  both  very  strong  energy  deposition  and  for  relatively  weak  energy  deposi¬ 
tions.  Strong  shocks  leave  large  yet  subsonic  velocities  in  their  wakes.  In  terms  of  the  ambient  sound 
speed  (at  the  pulse  location),  this  velocity  of  expansion  may  be  quite  supersonic  making  our  approxi¬ 
mation  actually  underestimate  the  vorticity  generated.  On  the  other  hand,  for  weaker  energy  deposi¬ 
tions  the  shock  speed  is  still  slightly  larger  than  cs  but  the  flow  speed,  and  hence  the  expansion  rate, 
can  be  much  smaller. 

These  weaker  cases  may  well  be  treated  as  dynamic  adiabatic  expansions  where  Ri/R^—  1  is 
small  but  detailed  numerical  simulations  will  be  needed  to  calibrate  the  improved  theoretical  model.  It 

may  well  be  that  the  effective  expansion  time  r  and  hence  the  effective  expansion  velocity  should 

involve  the  maximum  of  the  shock  transit  time  from  to  /?o  and  the  rarefaction  wave  transit  time  to 
the  axis  from  /?0.  The  effect  of  reducing  the  vorticity  generation  rate  at  low  deposition  power  means 
that  there  will  probably  be  a  shift  of  the  least  mixed  channel  to  larger  size  since  not  only  R{  -  R0  but 

also  the  effective  expansion  velocity  -y-  decrease  when  the  overpressure  T  is  small. 

In  reality  the  energy  deposition  will  probably  depend  on  the  existing  channel  profiles.  Lower  den¬ 
sity  often  means  reduced  energy  deposition  although  perhaps  equal  or  near  equal  local  heating  may 
result.  Furthermore,  circular  energy  deposition  will  no  longer  apply.  While  it  may  not  yet  be  possible 
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to  evaluate  the  channel  effect  on  subsequent  pulse  propagation,  the  effect  on  the  channel  of  an  energy 
deposition  with  given  cross  section  should  be  easier  to  determine.  It  will  be  particularly  important  to 
characterize  these  phenomena  quantitatively  because  a  reduction  in  the  beam-generated  turbulence  is 
likely.  This  is  turn  reduces  the  turbulence  level  and  hence  convective  channel  cooling.  For  near  term 
experiments  these  refinements  and  the  corresponding  reduced  cooling  for  low  intensity  beams  will 
probably  help  raise  the  average  channel  temperature.  For  more  intense  beams,  the  turbulence  problem 

Um 

may  well  be  worse  than  we  predict  here  because  -j-  may  be  larger  than  cs. 

Distortion  of  either  a  two-dimensional  or  a  three-dimensional  variety  will  also  contribute  to  the 
vorticity  and  hence  turbulence  generation.  These  additional  asymmetries  must  also  be  analyzed 
theoretically.  Even  in  the  initially  2D  cases  three-dimensional  effects  are  likely  to  appear  and  reach  full 
amplitude  by  the  time  an  averaged  description  of  the  "turbulent”  channel  cooling  will  become  useful. 
This  is  not  likely  to  occur  before  a  few  pulses  at  least  so  the  use  of  detailed  (and  probably  expensive) 
calculations  is  necessary  to  bridge  the  gap  between  a  single  pulse  interacting  with  a  previously  bored 
smooth  channel  and  a  roughly  circular  turbulent  cylinder  being  driven  by  a  repetitively  pulsed  energy 
source  near  its  axis. 

Questions  about  the  vorticity  generation  in  highly  distorted  geometries  must  be  resolved.  As  the 
symmetry  plane  jet  developes  in  the  simplest  case,  the  hot,  double  channel  becomes  bisected  into  two 
channels  with  embedded  vorticity.  A  third  pulse  might  be  expected  to  bisect  these  two  channels— if  the 
location  and  timing  of  the  pulse  were  correct.  It  is  clear  that  the  timing  is  important  because  the  vorti- 
city  should  add  essentially  linearly  when  the  pulses  are  separated  by  a  time  appreciably  shorter  than  the 
mixing  time  estimated  above.  By  separating  the  pulses  sufficiently  in  time,  on  the  other  hand,  interest¬ 
ing  nonlinear  effects  will  result.  Proper  phasing  and  location  of  the  pulses  can  certainly  cancel  existing 
vorticity  in  some  regions  and  enhance  it  in  others. 

As  the  number  of  pulses  increases  so  also  will  the  geometric  complexity  of  the  resulting  compo¬ 
site  channel.  A  turbulent  situation  is  expected  to  result  in  which  the  migrating,  interacting  vortex 
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filaments  shear  and  stir  the  hot  expansion  channels.  This  mixing  brings  cold  entrained  air  into  proxim¬ 
ity  with  the  hot  channel  air  on  spatial  scales  short  compared  to  an  expanded  channel  radius  and  thus 
molecular  thermal  conduction  cooling  is  expected  to  be  enhanced.  Questions  in  the  turbulence  regimes 
occur  at  a  number  of  points.  Where  the  beam-generated  vorticity  goes  and  how  it  evolves  are  clearly 
controlling  issues.  The  same  reduction  of  scales  by  the  expansion  induced  jets  which  tends  to  favor 
enhanced  molecular  thermal  conduction  will  also  favor  enchanced  viscous  dissipation  of  the  driving 
vortices  themselves.  Since  the  decay  of  vorticity  reduces  the  convective  cooling  rate,  a  balance  will  be 
obtained  eventually  where  the  production  of  vorticity  is  equal  to  the  dissipation  rate. 

In  standard  turbulence  phenomenologies  there  is  no  good  way  for  dealing  with  transient, 
multispace-scale  phenomena  so  a  model  for  these  effects  will  have  to  be  developed  which  also  allows 
the  complex  channel  to  be  cooler,  on  average,  at  the  edges  than  in  the  center.  The  importance  of  tem¬ 
perature  fluctuations  also  cannot  be  overemphasized.  Knowing  the  average  channel  temperature  is  not 
sufficient  if  one  very  hot  region  takes  on  unusual  importance.  Thus  a  transient  turbulence  model  will 
also  have  to  be  constructed  which  allows  molecular  cooling  so  the  temperature  fluctuations  about  the 
mean  can  be  analyzed. 

Work  is  currently  underway  to  extend  this  analysis  to  the  other  types  of  asymmetry  and  to  per¬ 
form  detailed  simulations  of  the  late  time  channel  cooling  to  calibrate  the  mixing  time  formula  as  well 
as  the  vorticity  generation  formula.  Three-dimensional  simulations  of  the  early  time  and  the  develop¬ 
ment  of  a  spectral  dynamics  turbulence  phenomenology  for  the  late  time  will  be  necessary  before  a 
truly  interactive  multipulse  scenario  can  be  simulated  where  the  beam  deposition-vorticity  generation- 
channel  distortion  feedback  loop  is  incorporated.  A  framework  for  the  turbulence  analysis  has  been 
developed,  and  will  serve  as  the  basis  for  our  improved  analysis  of  these  issues  in  future  research. 
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Table  I.  Values  of  /2  from  numerical  simulations.  The  displacement  of  the 
second  pulse  from  the  channel  center  is  varied  (c  a  x<js0)  and  two  different 
estimates  of  Um  are  compared.  The  sound  speed  in  the  original  channel  at  the 
second  pulse  center  is  indicated. 


Um  =  C,(X0) 

Um  (eq.46) 

-f  2  (eq.46) 

-V(r) 

8.2  x  104cm 
sec 

0.10 

8.0  x  104cm 
sec 

0.10 

1.9  x  104cm2 
sec 

5.2  x  io4 

0.18 

6.3  x  io4 

0.14 

2.1  x  io4 

1.0 

3.9  x  104 

0.23 

5.9  x  104 

0.15 

2.1  x  io4 

1.5 

3.6  x  104 

0.19 

4.2  x  104 

0.17 

1.6  x  104 

2.0 

3.5  x  104 

0.15 

3.2  x  104 

0.16 

1.2  x  104 

T/tam 
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Time  after  first  pulse  (msec) 

Fig.  1  —  Measure  of  temperature  rise  versus  time  for  multipulse  channel  formation  using  Astron  (by  permis¬ 
sion  of  the  authors).  The  'classical  theory”  curve  indicates  that  cooling  occurs  faster  than  can  be  accounted 
for  by  thermal  conduction  alone. 
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Fig.  2  —  Schematic  of  the  initial  conditions  in  which  a  second  beam  pulse  of  radius  R0  propagates  parallel  to  the 
hot  channel  of  characteristic  radius  S0  formed  by  a  previous  beam  pulse.  The  second  pulse,  offset  from  the  first 
pulse  by  X0 ,  expands  to  radius  /tt. 
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Fig.  3  —  Symmetric  distribution  of  localized  vorticity  appears  as  two  extended  vortex  filaments  V  separated  by  a  dis¬ 
tance  2 y.  Flow  velocity  induced  by  these  vortices  is  shown  along  the  x  axis  halfway  between  the  filaments.  The  mixing 
time  Tm„  is  defined  here  as  the  time  it  takes  a  fluid  element  at  -S0  to  reach  S0  and  hence  effectively  bisect  the  hot 
channel. 
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Fig.  4  —  Integrated  vorticity  form  factors  assuming  a  "super -Gaussian"  density  profile.  Larger 
variation  with  b  is  observed  in  the  upper  figure  than  with  a  in  the  lower  figure.  Peak  vorticity 
generation  always  seems  to  occur  in  the  vicinity  of  .Y0  -  S0  and  -f2  ranges  from  0.1  to  0.4  in 
maximum  value  with  0.3  being  a  good  working  value. 


Fig.  S  —  Integrated  vorticity  form  factors  assuming  a  Bennett  density  profile.  Larger  variation 
with  b  is  again  observed  than  with  a.  Note  that  appreciable  form  factors  at  Xq  ~  0  2  5n  can 
be  obtained  when  b  is  small. 
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Fig.  6  —  Comparison  of  integrated  vorticity  form  factors  for  three  distinct  density  profiles  at  a  -  b  —  I/n/1  Relatively  little 
variation  in  the  form  factors  is  seen  for  the  widely  differing  density  profiles 
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Fig.  7  —  Initial  conditions  for  numerical  simulations  of  beam-channel  vorticity  generation  mechanism.  Only  the  upper 
half  plane  is  simulated  with  outflow  boundary  conditions  at  X  —  0  cm.  X  —  30  cm  and  y  -  7  cm.  A  5:1  grid  stretching 
has  been  employed  for  the  outermost  zones  to  move  the  boundaries  far  from  the  active  region. 
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Fig.  8  —  Plot  of  8y  versus  cell  number  for  ihe  variably  spaced  y  grid.  Adaptive 
gridding  reduces  possible  boundary  reflection  of  shocks  to  a  small  error  by  put¬ 
ting  boundaries  of  the  computational  domain  far  from  the  hot  channel. 


CELL  NUMBER 


Fig.  9  —  Plot  of  tx  versus  cell  number  for  the  variably  spaced  x  grid.  Only  cells  51-100  are  shown  because  the  spacing  from 
cell  1  to  cell  SO  is  the  mirror  image  of  ihe  spacing  from  $1  to  100. 
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Fig.  10  —  Symmetry  plane  pressure  and  density  profiles  for  the  case  where  the  second  pulse  and  first  pulse  chan¬ 
nel  centers  are  separated  by  1.0  cm.  The  first  pulse  characteristic  channel  radius  is  S0  -  1.0  cm  and  the  second 
pulse  characteristic  channel  radius  is  0.4  cm  (see  Figs.  2  and  7).  Data  points  are  shown  near  the  channel  centers  to 
indicate  good  resolution  used  in  the  fine  zone  region. 
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S0  =  1.0  cm 


Fig.  1 1  —  Time  history  of  the  integrated  vorticity  in  the  upper  half  plane  for  our  standard  cose.  Overall  vorticity  is  zero,  of 
course,  because  an  equal  strength  but  oppositely  directed  vortex  filament  is  generated  in  the  lower  half  plane. 
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Fig.  12  —  Schematic  diagram  showing  how  initial  vorticity  transient  arises  physically.  When  the  outwardly  expanding 
fluid  eventually  decelerates  to  rest,  most  of  the  transient  vorticity  cancels  because  the  acceleration  pressure  gradient 
must  change  sign  to  stop  the  expansion. 
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Fi*.  14  —  Integrated  vorticity  P(t)  for  three  cases  with  X0  >  S0. 
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Fig.  16  —  Comparison  of  numerically  determined  form  factors  (•  and  A)  with  the  theoretical  calculation  using  the  vorticity 
generation  formulation  and  setting  a  -  0.4,  b  -  1.0.  Agreement  is  good  in  both  magnitude  and  sign. 
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Fig.  17  —  Mixing  time  estimates  as  a  function  of  pulse  overpressure  I"  and  initial  beam  radius  R0.  Ut  “  103 
cm/ sec  is  generally  characteristic  of  air  at  STP.  The  formulae  summarized  in  Section  V  were  used  with  l/jl  "  y 
and  C,  «  5  x  10*  cm/ sec. 
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